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摘 要 : 为 了 研究 碳纤维 类 高 温 防 隔 热 复合 材料 中 热 辐 射 传输 特性 及 对 其 传 热 性 能 的 影响 ， 采 用 数值 模拟 的 方法 ， 针 对 2 维 
(2-D)，2.5 维 C.5-D) 编 织 的 碳纤维 复合 材料 的 光谱 辐射 特性 进行 了 研究 。 通 过 对 比 时 域 有 限 差分 法 (FEDTD) 和 Mie 氏 理论 计 
算 的 微米 级 尺度 下 单 根 无 限 长 圆柱 的 散射 强度 分 布 ， 验 证 FDTD 计算 碳纤维 热 辐射 特性 可 靠 性 的 基础 上 。 采 用 FDTD 分 析 了 
编织 方式 、 纤 维特 征 尺寸 、 纵 向 层 数 、 纤 维 间隔 及 基质 折射 率 对 其 辐射 性 质 的 影响 。 研 究 发 现 ， 均 质 块 材 与 实际 编织 结构 的 
块 材 光 谱 辐 射 特性 差异 明显 ， 散 射 强度 分 布 差异 可 达 2 个 量 级 。 
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The Research of Structural Influences on Carbon Fibrous Composite Material 


Radiative Spectral Properties 
YANG Liu XIE Ming AIQing 
(School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China) 


Abstract: In this work, aim to study carbon fibrous composite insulation materials” radiative spectral 
properties and the influences on heat transfer, using numerical simulation to calculate 2-dimension(2-D) and 
2.5-dimension(2.5-D) micro-structural influences on carbon fibrous composite material spectral attenuation 
properties. The reliability of the result obtained from FDTD is proven by comparing the results of the 
scattering intensity distribution of the single carbon fiber with the Mie Theory approximation. Based on 
FDTD, this paper analyzes the influence of diameters, longitudinal structural layers, interval in horizontal 
and vertical as well as the back ground refractive index on the spectral attenuation properties of carbon fiber 
braided structure both in 2-Dimension (2-D) and 2.5-Dimension (2.5-D). We also discuss scattered 
differences in stagger and aligned arrangements respectively. Finally ,it can be shown that different 
structures have their special radiative spectral properties, meanwhile, it is not accurate to regard the carbon 
fiber braided structure as the equivoluminal block, the maximum value of the different scattered intensity 
may up to two orders of magnitude. 
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0 引言 


碳纤维 复合 材料 一 种 性 能 优异 的 隔 热 材料 ， 常 用 ”受到 了 越 来 越 多 的 关注 ， 圆 柱 的 位 置 以 及 排列 的 形状 
于 高 温 热 防护 。 通 常 碳纤维 材料 由 微米 级 的 碳纤维 编 “都 会 影响 入 射 光 谱 的 散射 方向 与 能 量 ，Lee Siuchun 等 
织 而 成 ， 难 以 基于 几何 光学 进行 研究 。 目 前 ， 学 者 们 ”人 对 平行 放置 的 圆柱 提出 了 求解 随机 放置 在 吸收 介质 
通过 使 用 米 氏 理论 ， 变 换 波 动 方程 形式 等 方法 求解 结 的 圆柱 散射 性 质 的 方法 9 "1 但 上 述 研究 主要 针对 2 
构 简 单 的 无 限 长 圆柱 的 散射 性 质 "…， 对 无 限 长 圆柱 的 甘 结 构 ， 而 且 没 有 考虑 实际 材料 编织 等 结构 形式 的 影 
研究 也 更 加 细 化 ， 对 放置 在 吸收 性 介质 中 的 单 根 圆柱 问 。 本 文通 过 使 用 时 域 有 限 差 分 法 ， 直 接 对 麦克 斯 韦 
的 散射 性 质 进行 研究 9H， 并 提出 了 对 光线 入 射 单 根 多 方程 进行 离散 与 求解 ， 研 究 了 不 同 碳纤维 编织 结构 的 
层 介质 的 通用 计算 方法 "LI。 多 根 圆柱 光谱 性 质 的 研究 。 辐射 光谱 性 质 ， 分 析 了 结构 对 复合 材料 热 辐 射 传 
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输 的 影响 ， 为 高 温 下 复合 材料 设计 提供 了 参考 。 
1 数理 模型 
1.1 物理 模型 
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图 1 2 维 碳纤维 计算 模型 


Fig.1 2-D carbon fibrous computational model 


图 1 为 2 维 尺度 下 碳纤维 的 物理 计算 模型 , 其中， 


圆 形 区 域 为 碳纤维 ，M 表示 水 平 层 数 ，N 表示 纵向 层 
a 
射 。 

1.2 基本 方程 


针对 碳纤维 粒子 的 辐射 特性 问题 ， 本 节 介 绍 了 时 
有 限 差分 的 基本 方法 ， 选 用 了 FDTD 计算 当中 的 
PML 边界 层 作为 几何 边界 的 处 理 方法 5 。 
麦克 斯 韦 方程 是 一 组 用 来 表征 和 描述 自然 界 宏观 
电磁 现象 的 方程 ， 其 中 微分 形式 基本 可 以 表征 大 部 分 
电磁 学 基本 原理 ， 形 式 如 下 : 
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V:D=p (1) 

VyV:.B=0 (2) 

Re (3) 
ot 

Ds (4) 
ot 


在 研究 材料 光谱 特性 时 ， 散 射 因 子 与 散射 分 布 强 
度 分 布 可 以 很 好 地 表征 它 的 光谱 特性 ,本 文选 用 这 两 
个 参数 作为 参量 ， 研 究 不 同情 况 下 二 者 的 变化 规律 ， 
由 FDTD 算法 可 计算 出 节点 处 的 入 射 ， 散 射 的 波 印 廷 
矢量 : 


S ,=Re(E )xRe(H.,.) (5) 
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由 波 印 廷 矢量 辐射 强度 的 关系 ， 在 一 个 周期 7 中 ， 节 
点 处 的 平均 辐射 强度 : 


Re(E,, ) x Re(H..,) (6 ) 
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5 人 SCDdi (7) 


设 定 光 源 垂 直 于 圆柱 轴 向 入 射 ， 当 基质 介质 为 非 吸收 


性 材料 时 ， 投 射 到 圆柱 上 的 辐射 强度 不 变 ， 所 以 能 量 
为 : 
WW, Lie “ 4 (8) 
40 为 圆柱 的 投影 面积 ， 散 射 能 量 头 
WW =L 中 Tadc (9) 
C 为 结构 区 域 的 边 长 ，L 材料 长 度 ， 上 标 c 表示 当前 
节点 处 的 散射 强度 。 散 射 因子 的 计算 式 为 : 
W 
1 = (10) 
Osos 到 


1nc 


1.3 FDTD 中 完全 匹配 吸收 边界 (PML) 条 件 
在 计算 域 的 边界 设置 完全 匹配 吸收 边界 来 截断 计 


算 。 事 实 上 ，PML 是 一 个 有 厚度 的 介质 层 ， 介 质 层 的 
波 阻 抗 和 模拟 区 域 波 阻抗 完全 匹配 ， 即 : 
O oO 
i (11) 
€ LH 
2 (12) 
€ LH 


下 标 x，y，mx，my 分 别 表示 二 维 坐标 下 不 同方 向 的 
电导 率 和 导 磁 率 ， 则 PML 中 的 电场 分 量 可 表示 速 为 : 


六 二 藉 op jd lt 
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其 中 , G 可 参阅 文献 [14],co 是 光 在 真空 中 的 传播 速度 ， 
0 是 入 射 方向 和 交界 面 的 夹 角 ，PML 边界 层 不 影响 入 
射电 磁 波 的 吸收 效果 。 在 处 理 纤维 边界 时 ， 为 防止 产 
生 阶 梯 近 似 ， 使 用 了 共 形 网 格 技术 1。 


2 编织 碳纤维 辐射 特性 分 析 


为 了 验证 本 文采 用 FDTD 法 的 可 靠 性 ， 针 对 无 限 
长 圆柱 体 ， 对 比 了 FDTD 法 与 Mie 理论 的 计算 结果 ， 
如 图 2 所 示 ， 结 果 显 示 二 者 吻合 良好 。 
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scattered intensity distribution 
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图 2 Mie 理论 与 FDTD 计算 的 散射 强度 分 布 对 比 


Fig.2 Comparison of the scattered intensity distribution between 


Mie approximation theory and FDTD method 


采用 FDTD 法 ， 针 对 复合 材料 2-D 和 2.5-D 两 种 
常见 编织 方式 进行 数值 模拟 ， 通 过 计算 不 同 条 件 下 材 


料 的 散射 因子 等 参数 来 对 比 研 究 其 各 自 的 热 辐射 


能 。 
为 了 控制 计算 变量 ， 假 设 碳纤维 材料 的 电磁 
均匀 且 各 项 同性 , 查阅 相关 数据 "9, 得 到 波长 在 2.5 


性 


生 质 


hm 


至 10 hm 时 膨胀 石墨 的 复 折 射 率 ， 并 把 它 当 做 碳纤维 


的 电磁 参数 使 用 。 
2.1 2-D 编织 碳纤维 辐射 特性 研究 


选用 波段 在 2.5 hm 至 10 hm 的 平行 光 作为 光源 ， 
计算 模型 如 图 1， 设 碳纤维 的 截面 半径 为 1 pm， 且 规 


定 矶 纤维 圆柱 的 水 平 间 隔 (LZ,〉 和 纵向 间隔 (ZL;) 
同一 个 计算 模型 中 为 定 值 。 如 果 没 有 做 特殊 说 明 ， 


在 
则 


取 间 隔 为 0 hm， 基质 折射 率 为 1， 结构 的 体积 斥 寸 为 


24x24 hum2。 


wavelength/nm 


3 2-D 不 同 横向 间隔 散射 因子 随 波 长 的 变化 规 得 
Fig.3 2-D structure Scattering factors of different 
horizontal intervals 
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图 3 表示 横向 粒子 间隔 LK 对 2D 编织 结构 散射 因 
子 的 影响 ， 四 条 曲线 的 变化 趋势 基本 一 致 ， 全 部 随 着 


入 射 波长 的 增加 ， 散 射 因 子 逐 渐 减 小 ， 下 降 速度 逐 


减缓 ， 相 同 波长 下 ， 间 隔 越 大 时 ， 散 射 因 子 越 小 。 纵 


簿 渐 


向 粒子 间隔 工 对 2-D 编织 结构 散射 因子 的 影响 关系 与 


横向 间隔 相似 ， 散 射 的 变化 趋势 也 基本 一 致 。 不 同 的 
是 纵 问 间隔 在 5.5 hm 之 后 的 波段 内 , 曲线 出 现 了 交叉 ， 


表示 在 这 个 波段 范围 内 ， 粒 子 间 距 已 经 不 是 影响 散射 


因子 的 主要 因素 。 
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图 42-D 不 同 纵向 间隔 散射 因子 随 波 长 的 变化 规律 

Fig.4 2-D structure scattering factors of different 

longitudinal intervals 


基质 折射 率 对 2-D 编织 结构 散射 因子 的 影响 关系 


十 | 


在 


图 5 中 给 出 ， 从 图 中 可 以 得 到 ，4 条 曲线 全 部 随 着 


入 射 波长 的 增加 ， 散 射 因子 逐渐 减 小 。 其 中 ，m1.1 


对 应 的 散射 因子 变化 区 间 最 小 ，ns=1.3 对 应 的 散射 因 
子 变 化 最 大 ， 另 外 ， 随 着 基质 折射 率 的 增 大 ， 散 射 因 
子 先 减 小 后 增 大 。 基 质 折射 率 对 2-D 编织 结构 散射 因 


于 


加 
维 
围 


增 
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的 影响 作用 很 大 ， 特 别 是 在 高 基质 折射 率 中 。 


5 T T T Tr T T T T 


wavelength/nm 


图 5 2-D 不 同 基质 折射 率 散 射 因 子 随 波 长 的 变化 规律 
Fig.5 2-D structure scattering factors of different back ground 
refractive indexes 


在 原 有 结构 尺寸 下 ， 通 过 加 大 它 的 纵向 宽度 来 增 
碳纤维 层 数 。 分 别 计算 了 N=6、8、10、12 层 碳 纤 
的 散射 性 质 ， 发 现 除 了 N=12 对 应 曲线 在 短波 段 范 
内 会 出 现 短暂 增 大 。 从 整体 来 看 ， 随 着 纵向 层 数 的 


兽 加 ， 散 射 因子 逐渐 变 大 。 碳 纤维 半径 对 2D 编织 结 


散射 因子 的 影响 关系 可 由 变换 的 波长 体现 ， 本 质 上 
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来 说 ， 在 用 电磁 学 的 方法 计算 材料 的 辐射 特性 时 ， 需 
要 考虑 的 变量 是 纤维 的 特征 长 度 上 5J， 由 纤维 的 半径 和 
入 射 波长 共同 决定 ， 所 以 ， 改 变 入 射 波 长 的 范围 ， 可 
以 同样 当做 是 改变 碳纤维 半径 的 模拟 计算 。 
2.2 2-D 碳纤维 错 排 排 列 方式 对 散射 特性 的 影响 

不 同 碳纤维 排列 方式 也 会 改变 材料 的 散射 性 质 ， 
本 文 建立 了 错 排 的 计算 模型 。 模 型 中 纤维 的 纵向 层 数 
为 三 层 ， 第 一 层 与 第 三 层 纤维 7 根 ， 中 间 层 为 6 根 ， 
所 有 纤维 均匀 分 布 。 碳 纤维 的 截面 半径 设 为 3.5 hm， 
纵向 间隔 为 固定 7nm， 横 向 间隔 5 作为 变量 。 


wavelength/nm 


图 6 错 排 不 同 横向 间隔 散射 因子 随 波 长 的 变化 规律 
Fig.6 2-D Structure Scattering factors of different horizontal 
intervals in stagger arrangement 


可 以 看 出 ， 不 论 是 错 排 还 是 顺 排 ， 间 隔 的 大 小 是 
否 变化 ， 散 射 因子 都 会 随 着 波长 的 增 大 而 减 小 ， 但 
般 情 况 下 ， 相 同 的 间隔 和 波长 下 ， 顺 排 的 散射 因子 数 
值 要 大 于 插 排 ， 当 将 横向 间隔 继续 增 大 至 11 pm 以 上 
时 ， 散 射 因 子 曲线 几乎 重合 ， 可 以 认为 此 时 间隔 影响 
很 小 ， 排 列 结构 方式 成 为 主要 影响 因素 。 

2.3 2.5-D 编织 碳纤维 辐射 特性 研究 


图 8 2.5-D 编织 结构 模型 
Fig.7 2.5-D carbon fibrous computational model 


2.5-D 编织 结构 是 一 种 常见 的 纤维 结构 , 模型 如 图 
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7 所 示 , 2.5-D 编织 结构 中 的 碳纤维 可 以 分 为 纬 纱 和 经 
纱 ， 经 纱 缠绕 纬 纱 ， 处 于 弯曲 状态 ， 但 经 纱 的 横 截 面 
依然 看 做 圆 形 ， 纬 纱 和 经 纱 之 间 互 相 紧密 排 布 。 边 界 
条 件 以 及 光源 的 设 定 与 2-D 结构 完全 一 致 。 其 中 ， 碳 
纤维 半径 1 um， 横向 纬 纱 间 隔 亏 =2 nm， 纵向 纬 纱 间 
隔 了 一 0 um， 基质 折射 率 nj=1.0。 

与 2D 编织 结构 一 样 ， 横 向 纬 纱 间隔 是 影响 2.5D 
织 结构 孔 际 率 的 因素 之 一 ， 为 了 研究 横向 纬 纱 间隔 
对 碳纤维 红外 辐射 特性 的 影响 ， 分 别 设 定 了 不 同 横 向 
纬 纱 间隔 四 种 结构 的 2.5D 编织 碳纤维 , 横向 纬 纱 根 数 
为 3， 纵 向 纬 纱 层 数 为 2。 


或 


3 4 5 6 7 8 9 10 
wavelength/nm 
图 82.5-D 不 同 横向 纬 纱 间隔 散射 因子 随 波 长 的 变化 规律 
Fig.8 2.5-D structure scattering factors of transverse 
weft intervals 


从 图 8 可 以 看 出 ， 散 射 因子 曲线 的 变化 趋势 是 
随 着 波长 的 变化 而 快速 降低 ， 这 与 2-D 模型 的 结果 相 
似 ， 值 得 注意 的 是 在 2.5-D 结构 时 ， 散 射 因 子 随 波 长 
的 变化 幅度 明显 下 降 ， 所 以 ， 可 以 得 出 当 纤维 经 过 复 
杂 编 织 后 ， 其 本 身 的 辐射 光谱 特性 有 了 改善 ， 降 低 了 
不 同 波长 对 材料 的 影响 。 


= 


scattered intensity distribution 
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图 92.5-D 编织 结构 和 等 体积 块 材 的 散射 强度 分 布 
Fig.9 The scattered intensity of bulk and 2.5-D woven fiber 


将 实际 编织 结构 的 纤维 材料 与 均 质 块 材 的 散射 强 


度 分 布 进行 对 比 , 如 图 9, 在 整个 散射 角 范 围 内 , 2.5-D 


结构 纤维 的 散射 强度 分 布 与 块 材 明 显 不 同 ， 某 个 角度 
上 数值 差 在 2 个 数量 级 之 上 ， 而 且 经 过 编织 后 ， 各 个 
方向 上 的 辐射 强度 要 比 均 质 块 材 更 加 均匀 。 
3 结论 

本 文通 过 数值 模拟 分 析 了 两 种 编织 结构 (2-D， 
2.5-D) Re 普 辐 射 特性 ， 比 较 了 编织 方式 、 纤 
维特 征 尺 寸 、 纵 向 层 数 、 纤 维 间 隔 及 基质 折射 率 对 其 
局 射 性 质 的 影响 。 研究 发 现 ， 纤 维 顺 排 时 ， 在 固定 的 


几何 区 域内 增加 碳纤维 个 数 则 会 使 


高 ; 


射 效应 也 越 强 。 在 相同 的 间隔 和 波长 
因子 数值 要 大 于 错 排 ， 但 不 论 是 | 
波长 的 增 大 , 最 后 都 趋 于 
和 2. 
碳纤维 用 复杂 形式 编织 起 来 ， 


材料 的 散射 因子 升 
基质 的 折射 率 与 纤维 折射 率 的 数值 差距 越 大 ， 散 
， 顺 排 的 散射 
顺 排 还 是 错 排 ， 随 着 
个 定 值 。 通 过 比较 块 材 2-D 
5-D 编织 结构 的 散射 强度 分 布 ， 发 现 通 过 把 多 根 
能 够 很 好 的 提升 材料 的 


光谱 稳定 性 。 


Pe 


生 质 


散射 强 


得 到 


因此 ， 不 同 的 结构 的 碳纤维 材料 的 辐射 罚 射 光谱 
有 其 特殊 性 ， 在 计算 时 ， 如 果 简 化 做 块 材 处 理 ， 
度 分 布 之 间 的 差异 可 能 达到 2 个 数量 级 以 上 ， 
的 计算 数据 并 不 准确 。 
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